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Abstract 

We consider finite temperature 1-loop diagrams with hard loop momenta and an ar- 
bitrary number of external gauge fields when the external momenta are either soft, 
or near the light cone and nearly collinear with the loop momentum. We obtain a 
recursion relation for these diagrams which we translate into an equation for their gen- 
erating functional. By integrating out the soft fields while keeping two collinear ones we 
find an integral equation, originally due to Arnold, Moore, and Yaffe, which sums the 
bremsstrahlung and pair annihilation contribution to the thermal photon production 
rate. 
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1 Introduction 



Particle production in thermal and non-thermal systems plays an important role in rel- 
ativistic heavy ion physics and in cosmology. Typically there is some weakly interacting 
particle species with a low abundance in a hot medium, from which these particles are 
then produced. Examples are photons produced in a quark-gluon plasma, or dark mat- 
ter candidates such as gravitinos, axions, or axinos produced in an otherwise thermal 
universe and at reheating after inflation. 

The production of photons from a thermal quark-gluon plasma is remarkably com- 
plicated. At leading logarithmic order, that is, at leading order not just in the strong 
coupling as, but rather in a s log(l/as) only 2 — >■ 2 scattering processes contribute 
and the rate has been computed by Kapusta et al. [I] and by Baier et al. [2]. How- 
ever, already for the complete leading order contribution the production cannot simply 
be understood in terms of scattering processes involving only a handful of particles. 
Bremsstrahlung and pair annihilation involve multiple interactions via soft gluon ex- 
change. Such processes are not suppressed, despite the large number of interactions 
[31 HI- The emission occurs almost collinearly so that internal lines are nearly on-shell 
and compensate the suppression. Viewed in position space the radiated photon and its 
source overlap over large distances, and the interference of different interactions cannot 
be neglected. This leads to the so-called Landau-Pomeranchuk-Migdal (LPM) effect 
[U El El E] • The complete leading order photon production rate has been computed at 
leading order by Arnold, Moore and Yaffe by summing all relevant diagrams. 

In this paper we present a different approach. We consider hard particles with 
momenta of order T which propagate through a gauge field background. The gauge 
field momenta are either soft (k ~ gT) or are almost on-shell with virtuality k 2 ~ g 2 T 2 
and almost collinear with the loop momentum. Formally we integrate out the hard 
modes in this soft and collinear background, leaving us with an effective theory for 
soft and collinear modes. The effective theory is described by an equation which has 
a similar structure as the non-abelian Vlasov equations [8] describing Hard Thermal 
Loops [9]. Then we distinguish the soft and collinear fields and identify them with 
the soft gluon field and the photon field, respectively. In a second step we integrate 
out the soft gluons to obtain an effective theory for the collinear gauge fields only. % 
The effective theory takes the form of an integral equation which has been obtained 

3 This approach is similar to the one used in Ref. [TUj , where soft gauge fields were integrated out 
to obtain an effective theory for ultrasoft (k ~ g 2 T) fields. 
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previously in the calculation of the photon production rate [I]. 

This paper is organized as follows. In Sec. |2]we briefly recall the relation between 
thermal photon production rate and the finite temperature polarization tensor. The 
main part of the paper is contained in Sec. [3]with the calculation of the 1-loop diagrams 
in a soft or collinear gauge field background. The key steps of our calculation are the 
approximation Q and the partial fractioning ( fit)]) . These are used both to compute 
the 2-point function (Sec. 13. 2j) and then to obtain the recursion relation for the n-point 
functions in terms of (n — l)-point functions (Sec. 13. 3p . Then all n-point functions are 
put together in an effective action (Sec. 13. 4p . In Sec. @]we integrate out the gluon fields 
and in Sec. Owe discuss the generalization from scalar quarks to spin-1/2 quarks. We 
summarize and conclude in Sec. El Finally, in Appendix |A] we show that the connected 
pieces which we encounter in Sec. H] vanish. 

Notation and conventions: We use the metric with signature H . 



2 The thermal photon production rate 

The rate per unit volume at which real photons with momentum k and polarization 
vector are produced in a hot QCD gas can be written as 

(27r) 3 2A; ^ = e%e* u Il^ (k) (1) 

where k° = \k\ is the photon energy and 

U<(k) = J d 4 xe ik - x (J"(0)r(x)) (2) 

is a 2-point function of electromagnetic current operators J M . We have not included 
the electromagnetic coupling constant in J M , which therefore appears in Eq. flTj). The 
2-point function ([2]) can be written in terms of the retarded polarization tensor 

IL%{k) = 2f B (k°)ImIlZ(k) (3) 

where /b(-B) = l/(e E ^ T — 1) is the Bose-Einstein distribution at temperature T. 

Among the leading order contributions to the production rate ([T]) are QCD Compton 
and qq annihilation. These are 2 — > 2 scattering processes for which only a few diagrams 
contribute. 

We will only be concerned with bremsstrahlung and pair annihilation contributions. 
There a quark undergoes multiple scattering with soft momentum transfer q ~ gT and 
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at some point emits a bremsstrahlung photon or annihilates with an antiquark into a 
photon. 

Therefore we want to evaluate IP" with a quark interacting via soft gluons only. For 
simplicity we first consider scalar quarks, the generalization of our method to spin- 
1/2 quarks is straightforward and is described in Secj5j The first step we will take is 
to integrate out the scalar quark fields in a soft or collinear gauge field background. 
That means that we have to calculate diagrams with two external photon lines and an 
arbitrary number of soft external gluon lines. 

3 1-loop diagrams with soft or collinear external 
gauge fields 

We consider thermal 1-loop diagrams with an arbitrary number of external gauge field 
legs. For the bremsstrahlung and pair annihilation contribution to photon production, 
two of the gauge field momenta are hard (k ~ T) and correspond to the produced 
photon. The remaining ones are soft (k ~ gT) gluons. The particle in the loop 
corresponds to a quark which suffers soft scattering via gluon exchange and radiates 
the photon. 

Even though the photons and the gluons play very different roles, our approach 
allows for a unified treatment, and in this section we do not have to distinguish these 
two. First we review the relevant kinematics which has been extensively discussed 
in Ref. jl]. It allows us to simplify propagators and vertices. We obtain a recursion 
relation between a diagram with n external gauge field lines and the difference of two 
diagrams with n — 1 external lines. Then we consider the current induced by the gauge 
fields which is the first derivative of the generating functional of all n-point functions. 
The recursion relation turns into a generalized kinetic equation. It can be viewed as a 
generalized Vlasov equation which contains a convective term and a force term. 

3.1 Kinematics and power counting 

The hard momenta we are considering are all almost collinear. Up to higher orders they 
all point into the same direction which we denote by v with v 2 = 1. The 3-momentum 
components in the indirection are denoted by 

P\\ = P ' v (4) 
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We define the light-like vector v = (l,v). One has to account for three distinct mo- 
mentum scales. 

1. The emitting charged particle, which corresponds to the particle in the loop, and 
the emitted particle both have p» of order T, which is our hard scale. 

2. All 3-momenta perpendicular to v are soft, p± ~ gT. Furthermore, all momen- 
tum components of the gluons are soft, ~ gT. 

3. Finally, all 4-momenta k are 'collinear', v ■ k = (ko — k») ~ g 2 T. 

This includes the case that the emitted photon is off-shell by an amount k 2 ~ g 2 T 2 
which is relevant for dilepton production. 

All momenta in the loop have k 2 ~ g 2 T 2 . Therefore the propagators are sensitive 
to the so-called asymptotic mass m ~ gT, which is given by the real part of the 
thermal self-energy, with hard loop momentum, of a light-like hard particle [TTJ . It is 
thus generated by integrating out the gluons with hard momenta. This procedure does 
not yet yield a thermal width which is only generated by integrating out the gluons 
with soft momenta. For scalars and fermions in the representation r of the gauge 
group 

m 2 = \c 2 {r)g 2 T 2 (5) 

with the quadratic Casimir of the representation r. In a SU(A) gauge theory C^r) = 
(A 2 - 1)/(2A). 

3.2 2-point function 

In this section we explicitly compute the 2-point function. It turns out that all n-point 
functions with n > 2 can be obtained from (n — l)-point functions through a simple 
recursion relation. Furthermore, all kinematic approximations which are needed for 
the general case already appear for the 2-point function. 

We work in the imaginary time formalism, where the loop integral consists of a sum 
over imaginary Matsubara frequencies and an integral over 3-momenta. The external 
momenta must also be taken imaginary and can be continued to real values only after 

It is oftentimes referred to clS TTIqq . 

5 Therefore, in contrast to [4] we never include the (IR divergent) thermal width in the quark 
propagators, only the asymptotic mass. 
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the sum has been performed. The Matsubara sum is performed as usual by writing it 
as a contour integral 



p =in2wT J C 



dpo 
2-ki 



\ + /b(Po) 



h( Po ) (6) 



where the integration contour goes up on the right of the imaginary axis and goes down 
on the left of the imaginary axis. Then we close the contour around the poles of the 
propagators. After evaluating the integral ([SD, we can continue ko from a Matsubara 
frequency towards the real axis. The retarded polarization tensor is obtained by taking 
ko = Re(fco) + ie. 

Now we describe the two main steps of our calculation. In the imaginary time 
formalism the scalar propagator appearing in Fig. [1] is 

A(p) = ~ l 2 = J (7) 

p z — m z v ■ p v ■ p — p ± — m z 

with v — (1, —v) . The other propagator in Fig. Ufa) will only be parametrically large, 
of order g~ 2 T~ 2 , when v ■ p ~ v ■ k because v ■ k ~ fi^T. In this case po approximately 
equals p\\ , and v ■ p ~ 2py . Therefore we can approximate 

A(p) ~ J-D(p) (8) 
2 P|| 

with 

= wp- (pi + m 2 )/(2 P ||) ^ 

After this approximation the propagator, viewed as a function of the energy variable, 
has only one pole. □ For the following it turns out to be very convenient to partial 
fraction the product of the two propagators, 

D(p)D(p -k) — — X — [D(p) - D(p - k)} (10) 
e{K,p) 



6 Note that we have made the approximation (J5J before the Matsubara summation, that is, when 
Po is still purely imaginary and of order T, even though is valid only for real frequencies. This 
is justified, however, because the second pole in po of the propagator, which gets lost through the 
approximation (jSJ), would give v ■ p ~ — 2p\\ ~ T. With this value of v ■ p the second propagator 
A(k — p) would be of order T~ 2 . This is suppressed relative to the contribution we keep, for which 
we have A(fc — p) ~ g~ 2 T~ 2 . The same type of argument applies to the poles of A(fc — p). 
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where 

2(p[| - fc||) 2p,| 

denotes the difference of the poles of the two propagators which is of order g 2 T. 

Both vertices are proportional to (2p — k)^. We associate the factor (2p||) _1 of Eq. (jSJ) 
with the left vertex in Fig. Ufa) rather than with the propagator itself and define the 
resulting vertex factor 

V^p,p-k) = ^-(2p-kr (12) 

2pi| 

The different components of V have different orders of magnitude. The largest com- 
ponent is in the direction of v and is 0(1). The transverse components are of order 
g, and the w-component is 0(g 2 ). For computing the production rate of real photons 
one needs the transverse components. It is therefore not sufficient to take into account 
only the leading order piece. We also include the transverse components, and the 
approximation 

V"(p,p- k )~^ [i 2 P\\ ~ fc llK + (2p - k)j_] (13) 

is implicitly understood in the following. 
We write the polarization tensor as 

n£(*) = / j^Vfap-k)* (t a fr u (k,p)) (14) 
which turns out to be convenient for discussing the general n-point function, and define 

^(P||,fc||) = /b(P||) - /b(ph - An) (15) 
Our result for the 2-point function is then 

n*(fc,p) = — p h k \\)V„(p-k,p)t b (16) 

As we already stated below Eq. (JSJ, we do not include the gauge coupling in IL^. 

There is also the tadpole diagram containing the 4-point vertex which contributes 
to the 1-loop polarization tensor. It does not have the collinear enhancement since it 
does not depend on the external momentum. But its transverse components are of the 
same order as those of the diagram in Fig. [1] (a). However, due to the independence on 
the external momentum the tadpole diagram has no discontinuity and therefore does 
not contribute to the production rate. Therefore we do not consider it here. 
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Figure 1: 1-loop diagrams with soft or collinear external gauge field lines. Only the 
2-point function (a) needs to be calculated explicitly. The n-point functions (b) are 
related to the (n — l)-point functions by a recursion relation. All external momenta 
are outgoing. 4-momentum conservation implies k = X^=i kj ■ 

3.3 Recursion relation for n-point diagrams 

Now we consider the diagram with n > 2 external gauge field lines in Fig. QJb). We 
will find that it can be recursively related to diagrams with n — 1 external lines. To 
obtain this relation we pick out one vertex, the leftmost one in Fig. DJb) which carries 
momentum k. The remaining ones have incoming momenta kj with j = 1, . . . , n — 1. 
As for the 2-point function we can use the approximation (jSJ) for the propagators. 
After the partial fractioning ([TO]) this diagram is proportional to the difference of the 
diagrams in which either the vertex with k\ or the vertex with is omitted. 

Each vertex carries a generator t a of some gauge group. As before we do not include 
the gauge couplings in the vertices. From the calculation of the 2-point function we 
know that in the vertices we have p° ~ p\\ . In analogy with Eq. (1141) we write 

n (n) ;:^(^ ■ ■ ■ , = / 4^ v Mp- [t a ^ {n) ^Zz\(h, k n „ uP 

(17) 

We only have to consider the two propagators D(p) and D(p—k) which are connected 
to the left vertex in Fig. QJb). We apply the same approximations as for the 2-point 
function including the partial fractioning f fTOj) . The two terms in Eq. fflOl) are then 
proportional to diagrams in which either the propagator with momentum p or the one 
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with momentum p — k has been omitted from Fig. 0Jb). Thus, if one also leaves out 
the vertex factors connected to these propagators each of the two terms gives a (n — 1)- 
point function. For the second term in Eq. ffTUl) we obtain a contribution in which 
the propagator with momentum p has been omitted. Here we perform a shift in the 
summation variable, p° — > p° + k®. Then the remaining propagators are the same which 
appear in the (n — l)-point function with incoming momenta fe, • • • , k n -i, but with the 
loop 3-momentum p replaced by p — k\. Therefore we can write H^ n \ki, . . . , k n -i,p) 
in terms of the difference of lT n_1 ) with either fci or k n _i omitted, 

= - n^^zlih: • • • , K^pW^Ap -k, P -k + K^-i 
+ v^(p - fci.p^n^^^Cfe, . . . , Up- fci) (is) 

3.4 The induced current 

Now we attach a gauge field = t a W^ to each vertex in Fig. [T](b) except to the 
one with momentum k, and then sum over all n. The result can be interpreted as 
the current J M which is induced by the gauge field background. It can also be viewed 
as the first functional derivative of the effective action or generating functional of our 
diagrams. Working with the induced current is more convenient than with individual 
diagrams because one does not have to worry about summing over all permutations of 
external lines. Furthermore, it can be used for integrating out the soft gauge fields. As 
in Eqs. (1TI|) and (fT7j) we write 

j;(k) = J -0- V ,(p,p-k)tr (t a J(k,p)) (19) 
The "unintegrated" current J is given by 

oo n-1 , „ , 4 , v / n-1 \ 

j(k, P ) = E n [J ^f w ^) ^ 5 [ k - E k A 

xn^ri^.-.^i.p) (20) 
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Using Eqs. (TIB]) and ffTHj) one obtains a relation for J(k,p), 

e(k,p)J{k,p) =F{p h k\\)V{p-k )P )-W(k) 

d 4 q 



(2tt)< 



J(k — q, p)V(p — k,p — k + q) ■ W(q) 
-V( P -q,p)-W(q)J(k-q,p-q) 



4 Integrating out soft gluons 

The photon polarization tensor which enters the production rate ([[]) can be obtained 
from the diagrams in Fig. fT] by identifying two external lines with photons and the 
remaining ones with gluons. Connecting the gluon vertices with propagators and inte- 
grating over the gluon momenta will generate precisely the ladder diagrams studied in 
[3]. In addition, it generates quark self-energy insertions with soft gluon loops, which at 
LO are purely imaginary and correspond to a thermal width. These two contributions 
by themselves would be IR divergent, but their sum is IR convergent. 

In terms of the current (TT9]) it means that one of the background gauge fields is the 
photon field and all others are gluons. The gluon fields are integrated out, leaving 
only the photon. Then the current is just IL^A^, from which one can read off the 
polarization tensor IP^. 

We therefore distinguish between external photon and gluon fields and G M , and 
write 

W» = A" + (22) 

The photon carries k\\ of order T, while the gluon field has q\\ ~ gT. In order to 
compute the photon polarization tensor we consider one external photon field. In J 
we only need to consider terms zeroth and first order in A^ } J = J + J l5 and what we 
need to compute is J\. 



7 Fourier transformed with respect to k this relation has a similar structure as the non-abelian 
Vlasov [8] equations from which one obtains the generating functional of Hard Thermal Loops [9]. 
The linear part of the convective term v ■ k in the Vlasov equation got replaced by e(k,p). The force 
term corresponds to the first term on the RHS. Finally, the non-linear part of the covariant convective 
derivative is replaced by the integral on the RHS. In fact, if one could neglect k and q relative to p, 
then the integral would be proportional to the Fourier transform of the commutator [v ■ W, J] . Unlike 
the Vlasov equation our equation is non-local with two sources of non-locality. One is the terms with 
fcii in the denominator, and the other is due to the fact that the second J in the integral depends on 
p-q. 
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In the equation for J the function J 7 vanishes at leading order. Therefore Jq is 
suppressed compared to J\ and it can be neglected in the equation for J\. Thus the 
equation for J\ takes exactly the same form as Eq. (12 1|) for J, with W replaced by A in 
the inhomogeneous term, and with the two W's inside the integral replaced by G. We 
keep only the leading order piece of the gluon vertex factors, so that V(p—k,p—k+q) ~ 
V(p — q,p) ~ v and thus 

e{k,p)J 1 {k,p) = F(p h k\\)V(p -k,p)- A{k) 

~ J 70 [Ji(k ~ q,p)v ■ G(q) - v ■ G{q)J x {k - q,p h p ± - q ± j\ (23) 

Now we would like to integrate out the gluon field. We denote the resulting current 
by (Ji). In order to see how it works we write Eq. ( 123]) schematically as J\ ~ A + GJ X , 
leaving out all terms and all factors which are not relevant for the present discussion. 
Now we iterate it once to obtain J\ ~ A + G(A + GJ\). The term GA vanishes when 
we integrate out the gluons. Thus we can drop this term and write J\ ~ A + GGJ±. 
Integrating out the gluons then gives (J\) ~ A + (GGJ\). The two gluon fields can 
either be contracted with each other, or with the other gluon fields in Ji, that is, 
(GGJi) ~ (GG)(Ji) + (GGJi)connncctcd- In Appendix 1X1 we show that the connected 
part vanishes at leading order and can therefore be dropped. Thus, by integrating out 
the gluons one obtains a closed equation for (Ji) of the form (Ji) ~ A + (GG)(Ji). 

Now we can become more explicit. After iterating the integral equation once and 
integrating out the soft gluons, they have disappeared as external particles and appear 
only in terms of their propagator, 

(GZ(q)Gl(q')) = g 2 S ab A^(q)WS(q + q') (24) 

Since we are interested in the electromagnetic current we put t a = 1 in Eq. (|T9l) . Using 
trl = d(r) and t a t a = C2(r)l, we obtain 



e(k,p) tr(Ji(fc,p)) = d(r)F(p h k\\)V(p-k,p)-A(k) 



2 . d A q vWA^q) 



~ 2C ^ J (27r) 4 v ■ (fc - q) tr [{MKp))- {Mk, PhP± - q ± ))\ (25) 

The square bracket neither depends on q nor on q\\ , and the integrals over q and q\\ can 
be performed. Eq. (j25p is thus an integral equation which determines the transverse 
momentum dependence tr(Ji). Inside the integral we have approximated e(k — q,p') 
v ■ (k — q), i.e., we have neglected the terms containing transverse momenta and 
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thermal masses, even though they are of the same order as the term we kept. This 
is possible because the gluon propagator separately depends on q and q\\, which are 
both of order gT, and not on their difference v ■ q, which is of order g 2 T. Therefore 
the terms we have omitted only contribute to a higher order shift of the integration 
variable q\\. 

To perform the integrals over go and q\\ in Eq. ( |25|) we use again the imaginary time 
formalism. It means that we replace the integral (2tt)^ 1 J dq$ by a sum over Matsubara 
frequencies like in Eq. Furthermore, it means that ko in Eq. f[2"o"j) has to be an 
(imaginary) Matsubara frequency. Only after performing the sum over q one can 
analytically continue ko towards the real axis. We are interested in the production rate 
which is proportional to the imaginary part of the retarded propagator. Therefore we 
have to give ko a small imaginary part, i.e., ko = Ke(ko) + ie. In Appendix [X] we show 
how we perform the Matsubara sum. Using standard results for the HTL resummed 
propagators [12] and the sum rule of Ref. [13] one then obtains 



</0 



2ir v ■ {k — q) 



4" 



-T 



q1 



rriT 



(26) 



where we were able to neglect the dependence on v ■ k ~ g 2 T. That is, the only 
dependence on v ■ k enters through the imaginary part of ko- The Debye mass squared 
enters through the gluon propagators and is given by 



Til 



D 
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■{2N + N S + N 



f> 



(27) 



for a SU(iV) gauge theory with N s complex scalars and Nf Dirac fermions. 
For the production rate ([T|) we need the polarization tensor which we write as 



iWfc) 



d 3 p 



Vf,(p,p-k)Il u (k,p) 



{2i 



The reduced polarization tensor U u is related to J\ through 



tT(J 1 (k,p))^U w (k,p)A v (k) 
Therefore it satisfies the integral equation 



(29) 



e(k,p)Il u (k,p) = d(r)J r (pn,k\\)V l/ (p - k,p) 



+ iC 2 (r)g 2 T 



d 2 q± 
(2tt) 2 



n„(fc,p) - U. u (k,p\\,p_i - q ± ) 



(30) 
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We finally want to show that our integral equation f[3"Ul) can be reduced to the one in 
Ref. [I] for the production of real photons. We choose v in the direction of k so that 
k± = 0. Then we have v ■ k — 0, and thus e(k,p) = (p\ + m 2 )k\\/\2p\\{p\\ — k\\)]. Only 
the transverse components of n„ contribute since the polarization vectors in ([1]) are 
purely transverse. If we define a new function / via 

n ± (fc, P ) = 2 _^ " f(k, P ) (3i) 

we obtain the equation of the same form as in [I] 
2p ± = ie{k,p)f(k,p) 



g 2 C 2 (r)T 



q\ q\ + ml 



[f(k,p)-f(k,p h p ± - qi _)} (32) 



(2tt 

5 Spin- 1/2 quarks 

So far we have always been dealing with scalar quarks to avoid technical complications. 
There is no need to redo the entire calculation for spin- 1/2 quarks. The results above 
are still valid up to a few modifications which we now describe. 
We now have to deal with the resummed fermion propagator 

p-E{p) 

in the high temperature limit when the zero temperature mass can be neglected. In 
the plasma rest frame the self-energy S(p) takes the general form 



X(p) = a(p)p + b(p) 1 ° (34) 

Therefore chiral symmetry is not broken by thermal effects, and the left- and right- 
chiral fermions propagate independently. Since ~ g 2 T, we have a(p) ~ g 2 . Thus 
a(p)p is small compared to the tree level contribution and may be neglected. For p ~ T, 
p 2 ~ g 2 T 2 the propagator can then be written as 

m - - --r^-f < 35 > 

p z — 2o(p)p u p z — m z 
Here we have neglected terms of order g 2 T in the numerator and we have identified 

2 

b (P) = (36) 
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at p 2 = m 2 , where m is the asymptotic thermal mass in Eq. (j5J) . Since the left and right- 
handed fermions propagate independently, one may consider the photon production 
from left-handed quarks only. The complete rate is then twice as large. Thus one can 
deal with Weyl instead of Dirac spinors, the vertices contain a^ 1 instead of 7 M , and the 
propagator fl35|) contains a -p instead of f. Up to terms of order g 2 T, which we neglect 
in the numerator, p^ is light-like. Therefore we can write 

cr • p ~ 2p\\r)(p)rf (p) (37) 

where r)(p) is a normalized eigenvector of p ■ <y jp\\ with negative eigenvalue. Thus we 
find, similarly to Ref. [I], 

S L {p) ~ vip)v i {p)D(p) (38) 

with the same D(p) as in Eq. OH]). Note that, unlike in Ref. [I], our result is valid for 
both signs of p\\. We associate the spinors rj(p) and rf(p) with the vertices on either 
side of the propagator rather than with the propagator itself. Therefore the vertex 
factor now reads 

V(p, p-k)=rf(p- k)a"r)(p) (39) 

instead of Eq. (fl2j) . For real photon production one needs V up to order g, 

V{p,p- k) =v» + V£ + 0(g 2 ) (40) 

In the helicity basis the transverse components V±_ are particularly simple. We choose 
v as the 3-direction. Then for V ± = {V 1 ±iV 2 )/V2 one finds 

V + = ^- + 0(g 2 ), V- = P - + 0{g 2 ) (41) 

p\\ - h p\\ 

The other difference compared to scalar quarks is due to the Fermi-Dirac statistics, the 
Bose-Einstein functions get replaced by Fermi-Dirac distributions, so that J 7 becomes 

^(P||,A;,|)=/ P (p||)-/ P (p||-A;||) (42) 

6 Summary and Conclusions 

In this paper we have obtained an integral equation (Eq. ()2ip ) which sums all thermal 
1-loop diagrams with an arbitrary number of soft or collinear external gauge fields. We 
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have applied it to compute the rate for real photon production by bremsstrahlung and 
pair annihilation in a hot QCD plasma. 

Compared to the original calculation of the photon production rate in Ref. [I] our 
approach is significantly simplified by the fact the calculation is done in two steps. In 
the first step we have integrated out the hard momentum modes at one loop. The 
resulting effective theory is summarized by Eq. fj2Tj) . It has a similar structure as the 
non-abelian Vlasov equation which describes the Hard Thermal Loops for soft external 
gauge fields. In a second step we have integrated out the soft gauge fields corresponding 
to gluons. This results in the integral equation obtained earlier in Ref. [I] which sums 
all leading order ladder and self energy contributions to the photon polarization tensor 
for hard on-shell photons, and which thus describes the Landau-Pomeranchuk-Migdal 
effect on thermal photon production. 

Our approach should easily allow for generalizations. The method can be adopted 
to the production of other particles than photons, e.g. the production of spin- 1/2- 
fermions. We also hope for a possible generalization to non-equilibrium situations, as 
they occur e.g. in heavy ion collision. 

Acknowledgments This work was supported in part through the DFG funded 
Graduate School GRK 881. 

A Connected pieces 

Here we show that only the disconnected parts contribute in the calculation of the 
photon production rate at leading order, as claimed in Sec. HI To simplify the discussion 
we can drop the dependence on spatial momenta and on m 2 , and replace e(k,p) by 
v ■ k. Furthermore we leave out the vector indices of A and G. All these simplifications 
do not affect our argument. 

Therefore, for the present discussion, we may study, instead of the full integral 
equation (|23|) . the simplified version 




(A.l) 



where we have introduced the compact notation 




(A.2) 
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One can easily write down the solution to Eq. (lA.ip 

Mk) = ^ £s (1 cw »-(*-k*» ) A (* - 6 *) (a - 3) 

In the connected part (G , GJi) connnccte d, which was dropped in Sec. HJ the gluon field 
G(k\) is contracted with some G(A;m) with M > 3. Thus (GGJi) con nnectcd contains the 
Matsubara sum 

11 1 

T ^ A(fcl) ^ . (A; - A:0 ^ • (A: - ^ - A; 2 ) ' ' ' v ■ (k - k x ■ ■ ■ - k M -i) (A ' 4) 

We use the spectral representation of the propagators 

The thermal sum can now easily be performed using (Q. Out of the M poles, only the 
one at k® = u which gives /b(^) — T/uj = 0(1/ g) contributes at leading order because 
the gluons are soft. At all other poles the Bose distribution function would be 0(1) 
and the corresponding contributions can be neglected. 

After performing all Matsubara sums, we may analytically continue k° towards the 
real axis and replace it by k° + ie where k° is now real. Then (IA.4I) turns into 

L Disc A(fci) 

27T2 v ■ (k — ki) + ie v ■ (k — k\ — fc 2 ) + is 

(A.6) 



v ■ (k — k\ k M -i) + ie 

Here we have to integrate over the region where v ■ k± is of order g 2 T. That means that 
ki is equal to ki\\ up to terms of order g 2 T and that we may replace 

Disc A(ki) t^- Disc A(h lh ki) (A.7) 
k" Ki|| 

without changing the leading order result. Now we are done because all poles in the 
integrand lie above the real k®-axis. Therefore we can close the integration contour at 
—zoo and we obtain zero. This proves that at leading order only the disconnected part 
contributes. 

We finally remark that if we contract G(k\) with Gfaz), which corresponds to the 
disconnected contribution (and to M = 2 in the calculation above), this argument fails 
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since the integrand does not fall off rapidly enough at infinity and the integral over a 
closed loop would give a non- vanishing contribution. In fact, in that case we obtain 
the result (1261). 

In Ref. [1] the gluons where integrated out by (i) including the width in the quark 
propagators, and (ii) using Feynman diagrams for the remaining contributions. There 
it was shown that the leading order contributions are due to ladder diagrams with 
uncrossed rungs, and that ladder diagrams with crossed rungs or vertex corrections 
vanish at leading order. Even though we have not checked it explicitly, it seems to be 
clear that these diagrams are part of the connected pieces. 
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